BCS ansatz, Bogoliubov approach to superconductivity and 
Richardson-Gaudin exact wave function 



M. Combescot(i), W. V. Po eosov'^^'^\ and O. Betbeder-Matibet'^-'^) 
'^^^ Institut des NanoSciences de Paris, 
Universite Pierre et Marie Curie, CNRS, 
4 Place Jussieu, 75252 Paris Cedex 05 
'^'^'^ Institute for Theoretical and Applied Electrodynamics, 

Russian Academy of Sciences, Izhorskaya 13, 125412 Moscow, Russia 
(Dated: November 22, 2011) 

A Bogoliubov approach to superconductivity shows that the wave function ansatz proposed by 
Bardeen, Cooper, SchriefTer just corresponds to the ground state of the Bogohubov Hamiltonian, its 
energy being the extensive part of the BCS Hamiltonian energy. Extension to an arbitrary number 
of pairs in the energy layer where the potential acts, moreover shows that, although this wave 
function ansatz — with all pairs condensed into the same state — definitely differs from the exact 
wave function obtained by Richardson and Gaudin, it still provides the correct ground state energy 
for arbitrary pair number and potential strength, as obtained by analytically solving Richardson- 
Gaudin equations. This extension also reveals the existence of a super dilute and a super dense 
regime of pairs, with a gap different from the usual gap, these regimes being of importance to fully 
control the density-induced BEC-BCS cross-over. 

PACS numbers: 



I. INTRODUCTION 

This year, we celebrate the centenary of the dis- 
covery, by Kammerling Onnes [1], of superconduc- 
tivity which is one of the most fascinating phenom- 
ena of Solid State physics: indeed, it is so at odd 
from usual understanding that it took half a cen- 
tury to reach some acceptable microscopic picture 
of the phenomenon. The first clue was to under- 
stand that, in spite of their charge repulsion, two 
electrons can attract each other via the ion motion 
[2]. Although very small, this attraction, as shown 
by Cooper [3], can produce a two-electron bound 
state when acting in a region where the density 
of states is finite. The next step was to note that 
fermion pairs being boson-like particles, these two- 
electron bound states can condense into a collec- 
tive state quite different from the "normal" elec- 
tron gas. Although SchriefFer kept claiming that 
this condensation is quite different from a Bose- 
Einstein condensation of elementary bosons [4] , the 
wave function ansatz used in the Bardeen-Cooper- 
Schricffer (BCS) theory of superconductivity [5] is 
just the grand canonical ensemble version of all 
pairs being in the same state — as seen more in 
detail below. To support Schrieffer's claim, it is 
however clear that a state reading as the product of 
two linear combinations of free electron pairs is dif- 
ferent from two individual products because, due 
to the Pauli exclusion principle between fermions, 
one free fermion pair state is "missing in the sec- 
ond pair" . This Pauli blocking effect increasing 
with pair number, the state reading as a product 
of N identical linear combinations of free electron 
pairs is for sure very different from individual 
products, as for N elementary bosons in a Bose- 



Einstcin condensate. 

Nevertheless, even if more and more free pair 
states are missing when the number of correlated 
pair states increases — through what we called 
"moth-eaten effect" [6] — it is still of importance 
to know if the picture of superconductivity re- 
sulting from a product of identical electron pairs 
still is an acceptable picture of the phenomenon. 
This picture comes from an ansatz which has, as 
formidable support, the fact that it leads to easy 
calculations with results in agreement with exper- 
iments. Nevertheless, it is well known that wave 
functions are quite tricky, different ones possibly 
giving same values for quantities as averaged as 
the energy. 

The fact that Cooper pairs are not elementary 
but composite bosons pushes us to question the 
widely accepted idea that they all are in the same 
linear combination of free-pair states, as for an el- 
ementary boson condensate, even if, clearly, some 
collective effect takes place in this phenomenon 
- as seen from the fact that the linear combina- 
tion in the dense BCS limit is definitely different 
from the one of a single pair found by Cooper [3]. 
This question becomes even more relevant when 
considering the work done by Richardson [7] and 
by Gaudin[8] a few years after the BCS milestone 
paper on superconductivity [5]. Indeed, they suc- 
ceeded to write the exact form of the A'^ electron 
pair eigenstates for the Hamiltonian considered by 
Bardeen, Cooper, Schrieffcr[5] — this Hamiltonian 
actually being one of the very few exactly solvable 
model. This exact A^-pair wave function reads as 

B\R,) ■ ■ ■ B^RnM , (1) 
where B'^{Ri) is the following linear combination 
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of up and down spin electrons, 

^'(^') = E^;f^4t«U,, (2) 

k 

Wk being equal to 1 in the layer where the potential 
acts — called later on "potential layer". The i?i's 
are solution of N coupled equations 

In view of the second term of the above equation 
which can be traced back to the Pauli exclusion 
principle [9], the i?i's must be different. As a re- 
sult, the exact wave function (1) differs from 
pairs being in the same linear combination, as in 
the BCS ansatz. 

Over the last half century [10], these Richarson- 
Gaudin equations stayed analytically unsolved for 
arbitrary pair number N and potential V. In or- 
der to better understand the effect of Pauli block- 
ing between composite bosons — which forces us 
to consider a fixed number of pairs feeling the po- 
tential and thus amounts to stay in the canonical 
ensemble — we tackled these Richardson-Gaudin 
equations again recently[6, 11, 12]. In contrast to 
traditional BCS theory which corresponds to fill 
half the potential layer, we considered an arbitrary 
filling of this layer. We first succeeded to get the 
i?i's solution of these Richardson's equations in the 
dilute limit of pairs, i.e., for N Nc where Nc is 
the pair number above which single Cooper pairs 
would start to overlap. Later on, wc found a way 
[13] to reach the A^-pair ground state energy E^^, 
i.e., the sum '^^Ri, without calculating the i?i's 
individually. En turns out to have a nicely com- 
pact A^- dependence, 

p I — a 

within undcrcxtcnsive terms in (N/p)^, where p 
is the density of states taken as constant in the 
potential layer, a = exp(— 2/pF) while Ei = 
2eFo— 2ficr/ (1— cr) is the single pair energy found by 
Cooper, ei?Q being the Fermi energy of the frozen 
core electrons, i.e., electrons which do not feel the 
potential, and VL the potential layer extension, of 
the order of twice a phonon energy. 

It has been a surprise to note that this energy 
fully agrees with the ground state energy obtained 
by Bardccn, Cooper, Schricffcr, through a mini- 
mization of the Hamiltonian mean value calculated 
with their wave function ansatz [5] — which funda- 
mentally takes the A^ pair states condensed into 
the same linear combination, while the exact wave 
function (1) corresponds to different linear combi- 
nations. Indeed, the BCS ground state energy is 



known to read in the weak coupling limit, cr ~ 0, 

as 

Ebcs^\p^\ (5) 

where, in this limit, the gap scales as A ~ VL^/a. 
This result agrees with Eq.(4) for half filling, i.e., 
A^ = p^/2, which is the pair number correspond- 
ing to a potential extending symmetrically over a 
phonon energy scale, on both sides of the normal 
electron Fermi sea. Even if cases of wave functions 
very different from the exact one while giving the 
same energy, are known to exist, the fact that the 
exact A^-pair energy does support the BCS wave 
function ansatz still is rather puzzling because in 
this wave function, the A^ pairs are condensed into 
the same state while in the exact wave function, 
the A^ pairs all are in different states as a result of 
Pauli blocking. 

Actually, there is another approach to super- 
conductivity which allows reaching some micro- 
scopic understanding of this ansatz, although not 
yet helping to make link with Richardson-Gaudin 
exact wave function. It is based on the Bogoliubov 
transformation. This transformation allows an ex- 
act diagonalization of a part of the original BCS 
Hamiltonian for this Hamiltonian extended to the 
grand canonical ensemble. The ground state of 
this partial Hamiltonian just corresponds to the 
BCS ansatz. We arc then left with showing that 
the remaining part of the BCS Hamiltonian brings 
a negligible contribution in the thermodynamical 
limit, i.e., when the change from canonical to grand 
canonical ensemble is expected to be small. In 
fact, such conclusions have already been reached 
by Boboliubov [14] through a splitting the whole 
Hamiltonian in two terms, a "mean-field" Hamil- 
tonian and a perturbation. This splitting is similar 
to the one wc here use. Bogoliubov then uses this 
splitting to estimate from below and from above 
the ground state energy of the system. Instead, 
we here use a perturbative approach to, in a di- 
rect way, show that the perturbative part of the 
Hamiltonian produce undcrcxtcnsive terms only, 
provided that the arbitrary scalars introduced in 
the Bogoliubov procedure are chosen properly. In 
Rcf.[15], Bardeen and Rickayzcn have used single- 
particle Green's functions, in a form given by Gal- 
itskii [16], to reach similar conclusions. We also 
wish to mention a paper by Mattis and Lieb [17], 
in which it is shown that appropriate boundary and 
continuity conditions on the exact wave function in 
the large-sample limit, can lead to the equations of 
BCS theory obtained within the BCS ansatz. 

The paper is organized as follows. 

In section II, we briefly recall the BCS model 
Hamiltonian for superconductivity and the wave 
function ansatz proposed by Bardeen, Cooper, 
Schricffcr with cmphazis on its A^-pair form. 
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In section III, we propose a very simple presenta- 
tion of the Bogoliubov approach to the BCS prob- 
lem which leads to a splitting of the BCS Hamil- 
tonian as Hbcs = Hb -\- W. At this stage, both 
Hb and W are not fully fixed, so that they depend 
on a set of arbitrary scalars. We then show how 
the "Bogoliubov Hamiltonian" Hb can be diago- 
nalized. This procedure enables us to find the Hb 
energy spectrum. 

In section IV, we discuss some properties of the 
"Bogoliubov Hamiltonian" eigenstates. In partic- 
ular, we relate its ground state to the BCS wave 
function ansatz. 

In section V, we construct a perturbation theory 
for the total Hamiltonian Hbcs by treating W as 
a small perturbation. Using this approach, we de- 
termine the optimal set of scalars introduced in the 
Bogoliubov procedure, which makes the W contri- 
bution undcrcxtensive. The Hb ground state then 
reduces to the BCS ansatz when these underex- 
tensive terms are neglected. This gives a strong 
mathematical support to this ansatz in spite of its 
major difference with the exact form of the BCS 
ground state obtained by Richardson and Gaudin. 

In section VI, we concentrate on the ground 
state of the Bogoliubov Hamiltonian Hb ■ We con- 
sider not only the standard BCS half-filling config- 
uration, but also arbitrary fillings of the potential 
layer. Although these arbitrary fillings are difficult 
to achieve experimentally, their analysis enables us 
to reach a deeper understanding of the physics of 
Cooper-paired states. They also allow us to estab- 
lish a close link between two classical problems, 
namely. Cooper problem and the BCS model of 
superconductivity. 

In Section VII, we study the Hb excited states 
and derive the energy gap for an arbitrary filling of 
the potential layer. We also show the existence of 
two regimes, an extremely dilute and an extremely 
dense regimes of pairs which have a gap different 
from the usual gap A. These two regimes should 
have to be considered for a complete understanding 
of the BEC-BCS cross-over along the line proposed 
by Eagles [18] and by Leggett[19]. 

In Section VIII, we present some concluding 
comments. We, in particular, come back to the 
major physical discrepancy which still exists be- 
tween the BCS wave function ansatz and the 
Richarson-Gaudin exact wave function. 



II. BCS WAVE FUNCTION ANSATZ 

We consider a system made of N pairs of differ- 
ent fermions a and /?. In the case of BCS super- 
conductivity, these fermions are up and down spin 
electrons. For fermions with same mass, the free 
part of their Hamiltonian reads in terms of their 



creation operators and h]^ as 

Ho = Y^ ek(aj,ak + 6|,&k) • (6) 

k 

In usual BCS pairing, fermion-fermion attraction 
is reduced to processes between zero-momentum 
pairs. So, the total Hamiltonian reads as Hbcs = 
^0 + Vbcs with 

Vbcs = - ^ ^k'kak/&Lk/&-kak ■ (7) 

k',k 

In order to end with an analytically solvable prob- 
lem, the Vk'k scattering will be ultimately taken 
in a separable form l^cJk"^k- However, let us first 
keep its general form Vk'k- 

We want to determine the ground state of these 
N pairs. In order to have the lowest possible en- 
ergy, they must enjoy the attractive potential as 
much as possible: so, these pairs must appear as 
products of zero-momentum free pairs with cre- 
ation operator 

El = aih\ . (8) 

We thus expect the ground state for N pairs to 
read as 

\i^f)= E ^^(ki,...,kA,)<...i?t^|o) , 

ki,...,kiv 

(9) 

with (ki, • • • , k^v) all different due to the Pauli ex- 
clusion principle. 

The handling of this Pauli blocking effect exactly 
is definitely difficult. A smart way to overcome this 
difficulty is to turn to the grand canonical ensem- 
ble, with a pair number not fixed, as proposed by 
Bardeen, Cooper and Schrieffer through an ansatz 
constructed on the idea that, electron pairs being 
boson-like particles, they are likely to condense all 
into the same state 

B^=Y,^^Bi. (10) 

k 

As a result, Bardeen, Cooper and Schrieffer pro- 
posed a wave function ansatz reading as 

\^BCs) = E I^^^'^IO) • (11) 

By noting that \(I)bcs) can also be written as 
\^BCs)=e'''\Q)^\{e^-<\0) 

k 

=n(i+^kSi)io), (12) 

k 

smce bI^\Q) =^ due to Pauli blocking, we end 
for lySk written as Vk/^k, and \4>bcs) taken as 
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rjlipBCs), where i] = l/Hk^k, with a wave func- 
tion ansatz \4'bcs) appearing in its usual form, 
namely, 



(13) 



Wc then enforce lukl 



1 in order to 



have {4'Bcs\4'Bcs) = 1- The (Mk,Wk) coefficients 
are ultimately determined by a variational proce- 
dure, through the minimization of the Hamiltonian 
mean-value {ipBcs\H\'ipBCs), as found in all text- 
books [sec for example Ref. [20]]. 



III. BOGOLIUBOV PROCEDURE 
A. Rewriting the BCS Hamiltonian 

It is actually possible to give a precise mathe- 
matical meaning to the wave function ansatz pro- 
posed by Bardeen, Cooper and Schrieffer, by using 
a Bogoliubov approach to superconductivity. A 
very simple way to present this approach is to in- 
troduce a set of yet arbitrary complex scalars Zk 
and to split the BCS potential given in Eq.(7) as 



Vbcs ^£ + Vb+W 



(14) 



The two-body character of the BCS interaction 
now appears through W given by 



W 



E^k'k(i3^-4')(i3k-^k) 



(15) 



k'.k 



The Zk scalars will be ultimately chosen to make 
the W contribution to the ground state energy 
small. 

By extracting W from Vbcs, wc are left with 
Vb being a one-body potential, which is the goal 
of the procedure. It reads 



--JI^k'B,^, +h.c. 



(16) 



where the prefactor Ak' depends on Vk'k and the 
arbitrary scalars Zk through 

Ak' = J2 ^k'k ^k . (17) 



The scalar £ in Eq.(14) is then given by 

£ = Y.Vk'kZ^,Zk . (18) 

k',k 

This scalar is real since Vk'k = ^kk' which follows 

from Vbcs = Vbcs- 

While Vbcs conserves the particle number, Vb 
and W do not conserve it separately. In order to 
possibly treat W as a perturbation independently 



from Vb , it is thus necessary to turn to the grand 
canonical ensemble. This leads us to introduce the 
Hamiltonian 



H = Hbcs - ^-^N , 



(19) 



where N is the number operator for fermions a and 

/3, 



^ = E(«k«k + folM . 



(20) 



In this grand canonical ensemble, the Hamiltonian 
H then splits as 



H = Hb+W 



(21) 



where Hb is the "Bogoliubov Hamiltonian" given 

by 

Hb=£ + Y. {^k(4«k + &^_k^-k) 

k 

-AkSt - A^Bk} , (22) 
in which we have set ^k = ^k ^ 



B. Diagonalization of the "BogoHubov 
Hamiltonian" 



The Hb Hamiltonian is a one-body operator. It 
is thus easy to diagonalize it by introducing new 
operators 



Ok = a;kaj^ + j/k^-k 
4fe-k + J/k«k 



&Lk 



(23) 



In the Ak = limit, i.e., in the absence of poten- 
tial, the {x, y) prefactors must reduce to Xk = 1 = 
x'y. and yk — = y{^. If we enforce these new oper- 
ators to have the same anticommutation relations 
as the original operators aj^. and b^_y^, we then find 
that (a;k,2/k) and {x[^,y]J must be such that 



while 



dk, aj^. 



\yL\' 



^k' ' 



= = ?/k x'y. + Xk 2/k 
the other three anticommutators, 



(24) 
(25) 

(26) 



= 



«L4 



ak,6! 



(27) 

being automatically fulfilled for the operators de- 
fined in Eqs.(23). 
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It will appear as convenient to rewrite Eq.(26) 



as 



(28) 



When inserted into Eqs. (24,25), this gives 1 



IXkr = ■ (29) 

As a result, we also have |ykP = Ij/kP- 
If we now write aj^. as in Eq.(23), we get 



aldk = |a;kpaj^.ak + |ykP&-k&^k 



VkxlBy, + h.c. 



(30) 



and similarly for 6Lj^fo_k- By noting that &-k^Lk — 

1 — feLk^-k, the sum of these two products of Bo- 
goliubov operators thus reads 



4ak + &Lk^-k |ykP + |2/k 



' |2 



-(I 



a^kl^ - lyknaUk 



;,|^-|yk|^jfoLk6-k 
Vkxl - y[^x'^ ) Bk + h.c 



.(31) 



Since the prefactors of a^^au and b\_y^b-i^ are equal 
due to Eq.(29), it becomes easy to see that the 
hamiltonian Hb given in Eq.(22) also reads 



Hb=Sb + J2 ^k(4«k + 6Lk^-k) 



(32) 



provided that we take ^k = £'k(^|a;kP — |ykP 
£^k(2|.TkP - l) and A* = -Ek{ykxl - y{^x'* 



-2iJk|a;k| ik- As a result. 



kkP = l-|ykP = i(i + ^ 



Ai = -tkLBk + e 



(33) 



It then follows from i?k and ^k being both real, 
that |ikP = |Ak|V(^^k + ek)^ Since |tkP = 
lykP/kkP is also equal to {Ey, - Ck)/(i^k + ^k), 
we end with 



i?k = A/e^ + |AkP 



(34) 



the square root being taken with a positive sign in 
order for £'k to go to ^k in the absence of potential, 
i.e., when |Ak| 0. 

The scalar £b in the "Bogoliubov Hamiltonian" 
of Eq.(32) is then given, using Eq.(18), by 

Sb^Y. ^k'k 2k' 2k + ^(a - Ei,) . (35) 

k',k k 



Through these detailed calculations, we already 
see that the Bogoliubov procedure allows us to 
recover a few important results of the BCS the- 
ory, originally derived using a variational proce- 
dure based on the BCS wave function ansatz. 



IV. EIGENSTATES OF THE 
BOGOLIUBOV HAMILTONIAN 

A. Ground state 

Equation (32) shows that the Hb ground state 
energy reduces to £b while the excited states read 
as products of operators aj^. and ^Lj^ acting on the 
Hb ground state. 

This ground state is, by construction, such that 
Sk or &_k acting on it gives zero. It is then easy 
to relate this ground state to the \(I)bcs) ansatz 
defined in Eq.(12). Indeed, 



aplc^BCs) = [x*pap + Vphl^j (l + f/JpS^ ^ |0scs(p)) 

= {y; + ^P^;)blp\<l>BCs{p)) , (36) 

where \4>bcs{p)) contains all the k states of the 
I'pBCs) product except k = p. In the same way, 

'*\n^ I 



b-pl^BCs) = [y'p* - ^p<ja^|</>Bcs(p)) . (37) 

We easily see that these two states reduce to zero 
provided that ipp is chosen such that 



= 3^ = (38) 

Xp 

which always is possible due to Eq.(28). Since 
(0BCs|0BCs) - Ylmi + ¥'k5k)(l + ^kBl)\0) 

k 

=n(i+i^kp), (39) 

k 

while, due to Eq.(38), 1 -I- I^SkP = 1 + |yk/a;kP = 
1/l^kPj the normalized ground state of the Hb 
Hamiltonian thus reads, within an irrelevant phase 
factor, as 



|6) EE 



Bcs) = ni^ki(i-!^snio) 

= Y[{ui. + VkBlm , (40) 



in which we have set 



Uk = Fkl = Fkl 

I I J'k I / I J'k 
•^k -^k 



(41) 
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From this, we get the pair number in the 
Hb ground state through the mean value of the 
fermion number operator defined in Eq.(20) as 



7V = i(6|iV|6)=5^|«k| 

k 



Excited states 



(42) 



If we now turn to the lowest excited states of 
the Hb Hamiltonian, we find, using Eqs.(23), that 
they arc given by 



at|0) 



p"p 



^^pypjaplOp) 



-«^|0p) 



(43) 



where |0p) corresponds to the |0) ground state 
given in Eq.(40) without the p state in the k prod- 
uct. In the same way. 



&Lp|6) 



;^^'p|Op) 



(44) 



the excited states (43) and (44) having the same 
energy £b + Ep. 

The above equations lead us to choose Xp and 
x'p both real, which always is possible in view of 
Eqs.(23). The Hb excited states then take a com- 
pact form, 



am 



a^,|Op) 



^LplO) = bl 



plOp) 



(45) 



We note that the p state of the a fermion ap- 
pears through the pair operator — aj^^Lp in the 
ground state |0), while in the excited state aj,|0), 
this p state appears through only; so, the effect 
of the excitation operator is to break the 
pair by removing its b^_p component. Similarly, in 

^lp|0), the a|, component of the B^ pair operator 
is removed. 

Using Eqs.(28,33,41), we find 




while Wk, a priori complex, is related to Ak through 

* /* I* Ak 

■f^k — — Wk — — ^a;ktk = — , , 

(47) 

These again are standard textbook results derived 
through the minimization of the Hbcs Hamilto- 
nian mean- value calculated within the BCS ansatz 
(13). 



V. EIGENSTATES OF THE TOTAL BCS 
HAMILTONIAN 

In the above procedure, we have essentially ex- 
tracted the unpleasant two-body part of the origi- 
nal BCS Hamiltonian, through W. The Hb eigen- 
state spectrum is then straightforward to find. In 
order for these eigenstates to have some connexion 
with the original BCS problem, the "perturbation" 
W must only bring a small contribution. We do 
have some flexibility to reach this goal through the 
set of arbitrary scalars Zk introduced in rewriting 
the BCS potential as in Eq.(14). 



A. First-order correction 

A good choice for these Zk scalars clearly re- 
quires the first-order correction to the ground state 
energy (0|iy|6) to be small. Using Eqs.(23), we 
can rewrite the original BCS operators (flkjiiLk) 
in terms of the new Bogoliubov operators (aj^., ^Lk) 
as 



= Mkflk + f^k^-k, 



(48) 



Since fikl 
find 



= 5_k|0) by construction, we easily 



Sk|0) = (^Uk^-k - I'kakj (ukfik + I'k^^k 
= (u^v^-vlBl)\0) , 



(49) 



where Bj^ = Sj^^l^j^ creates a pair excitation. As a 
result, the first order correction to the Hb ground 
state energy £b induced by the W perturbation 
given in Eq.(15), reads as 



{0\W\0) 



k',k 



X (6|(uk'Wk' - -'k' 

X (ukfk-^k-i'k^k)|6). (50) 

Since (0|i?k'-S^|0) reduces to (5k', k, this gives 

(OIW^IO) =-^Vkk|«k|' 

k 

-"^Vi^'wiuk'Vk' - ^k'j (t^kWk - Zkj- (51) 



k',k 



By turning to the continuous limit, each sum over 
k brings a sample volume L^, while the potential 
matrix element Hc'k depends on sample volume as 
As a result, the double sum over (k',k) in 
(0|W^|0) is proportional to the sample volume; it 
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is thus extensive while the first sum is sample vol- 
ume free, so that it brings an underextensive con- 
tribution to the ground state energy £b which is 
negligible in the thermodynamical limit. In order 
to neglect the whole (0|W^|0) contribution to the 
ground state energy, we are thus led to cancel the 
second term of Eq.(51) by choosing the scalars Zk 
as 



2k = WkUk 



(52) 



As a result, Eqs. (46,47) give Ak', defined in 
Eq.(17), as 



Ak' = 

1 

1 ^ 



k^k = ^ Vk'kUk-yk 

k k 



, Vk'k- , 



(53) 



To go further and get some explicit results, we 
must choose a particular form for the Vk'k poten- 
tial. By taking it as separable, 

Vk'k = Vuji^'UJi^ , (54) 

with = Wk, the solution of Eq.(53) reduces to 

Ak = Acjk , (55) 

where A is solution of the so-called "gap equation" : 

^ Wk 



2 
V 



(56) 



The second order correction to £b is given by 



42) = (0|M/P^ 



1 



Pi_W\b) , (59) 



where Pi is the projector over the subspace per- 
pendicular to |0). Since the intermediate states are 
either B'^\0) or B^,B^\Q), this second order term 
splits as 



'^z^ ? — — ■ (.""i 



-Ep — Epi 



(2) 

The first term of £g has two V's and one p sum 
only so that it depends on sample volume as 
while the second term has two Vs but two p sums 
so that it is sample volume free. As a result, this 
second order term again gives a correction to the 
£b ground state energy smaller than L^, i.e., un- 
derextensive and thus negligible in the thermody- 
namical limit. 

By counting the number of potentials V and 
the number of p sums, it is possible to show that 
the higher order terms of the W expansion of the 
ground state energy give underextensive correc- 
tions. Consequently, within the choice of Zk scalars 
made in Eq.(52), the extensive part of the ground 
state energy for the BCS Hamiltonian in the grand 
canonical ensemble, indeed reduces to the scalar 
£b given in Eq.(35). 



B. Higher-order corrections 

Within the choice Eq.(52) for the Zk scalars, we 
find 

bI-z^ = uIbI - vt^B^ 

-UkWk(^a[ak + &^k^-k) • (57) 

Consequently, the W operator acting on the Bo- 
goliubov vacuum |0) gives rise to a state which 
contains zero, one and two pair excitations: 



21/^c^kWkWk|i'k|'^i|6) 

k 

vY,u^,u^ul,vlBlBl\0) . (58) 



k',k 



The above equation readily shows that the first 
order correction to the ground state energy £b^ 
namely, £^j^^ = (0|W^|0), reduces to the underex- 
tensive term — V^j. Wkl'^kl^, as already found in 
Eq.(51). 



VI. EXTENSIVE PART OF THE BCS 
GROUND STATE ENERGY 

In view of Eq.(19), the ground state energy 
of the BCS Hamiltonian Hbcs is related to £b 
through £b = £b + 2/iAf, where N is the mean 
value of the fermion-pair number in the system 
[see Eq.(42)]. The scalar Zk's arc chosen accord- 
ing to Eq.(52) for W to give an underextensive 
contribution to the ground state energy. We more- 
over take the potential in the separable form of 
Eq.(54) to possibly reach analytical results. Using 
Eqs. (17,55), wc then find 



Vk'k^k'Zk = V ^ Zk'Wk' E 2k Wk 

k',k k' k 



V 



(61) 



This gives the BCS ground state energy £b as 
V 



1^12 

£b=£b + 2fiN = ^ + 5I(^k - Ek)ujk + 2fiN 



(62) 

Note that a factor Wk can be introduced in the k 
sum because ^k — Ey = for Ak = 0, i.e. for 



8 



Wk = 0. This will allow us to split the k sum and 
calculate the two parts separatly. 

Due to Eqs.(42),(46) and (47), the chemical po- 
tential ^, hidden in ^k's, is related to the pair num- 
ber through 



A. Standard BCS configuration 

In the standard BCS configuration, the potential 
layer extends between epo and epg + ^ in a region 
where the density of states can be taken as equal 
to a constant p. Moreover, the potential is said 
to extend symmetrically with respect to the nor- 
mal electron Fermi sea which corresponds to pairs 
filling half the potential layer, i.e., TV = pQ/2. 

For half-filling and a constant density of states 
P; Eq.(63) gives, when turning to the continuous 
limit 



(64) 



It is straightforward to check that this equation is 
fulfilled for = ep,, + n/2. 

The gap equation (56), for a chemical potential 
set in the middle of the potential layer, then gives 



pd^ , § + Vt^ + a^ 

= plog 



This leads to 



A 



1 - e-2/pV 1 _ a ' 



(65) 
(66) 



with (T = e-'^/P^ . 

For ^ = epo + the two sums in Eq.(62) 

reduce respectively to ^k<^k = and to 



/n/2 
-0/2 



pn kT^ ~ A2 

as obtained through an integration by parts. If we 
now insert this result into Eq.(62), we end with 
a ground state energy in the canonical ensemble 
given by 



Since A = for 1^ = 0, the ground state energy 
in the absence of potential reduces to £b{V = 0) = 




2NeFo+ — P^'^ ^"^j the condensation energy 
induced by the attracting potential V is given by 



£BiV = 0)-£i 



(69) 

For A given by Eq.(66), ^/l + AA'^/Q'^ reduces to 
(1 -I- cr)/(l — (j); the condensation energy for half 
filling thus ends by reading as 

Uv ^0)-Sb = ^^ = |p a^(i - -) . 

(70) 

This result fully agrees with the expression of the 
A^-pair energy we have obtained [13] by analyti- 
cally solving Richardson-Gaudin equations in the 
canonical ensemble, namely, 

N(N- 1) 1 + cr 

En = NEi + , 71 

p 1 — cr 

the single pair energy obtained by Cooper being 
El = 2eFo— 2ricr/ (1 — cr). Indeed, the above equa- 
tion gives the condensation energy, i.e., the part of 
En which cancels when a = 0, as 

En{V = 0)-En^n(^1-^^ n . (72) 

We can check that this condensation energy re- 
duces to Eq.(70) for N = pil/2, whatever a, i.e., 
not only in the weak coupling limit, ct ~ 0. 



B. Arbitrary filling 

We now consider an arbitrary filling of the po- 
tential layer, namely, values of N not exactly equal 
to half the number of states between ep^ and 
epg + fl. In this case, we must solve two coupled 
equations (56,63) for the gap A and the chemical 
potential p. After switching from sum to integral, 
the integration limits being eFa~ P and ej^g — /Lt-t-fi, 
these two coupled equations read as 



2 _^^ £F„-M + » + V(£Fo-A^ + »)' + A^ 
PV eF„ -p+ v/(eFo - A*)' + A2 



and 



(73) 



2N 



^i^;~j7TWT^ - vi(-F„-pr+A^. (74) 

These equations can be solved analytically. Their 
solutions for A and fi read as 



pn\ pVLjl-a 



(75) 
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N f 2N 
p \ pil J 1 — a 



(76) 



It is easy to check that this solution reduces to the 
values of A and p obtained for half-filling, i.e., for 

pn/2. 

By inserting these A and p into the BCS ground 
state energy £b given in Eq.(62), we get 



£b = 2Nei 



N 



-2Nn 1 



N 



pQ J 1 — a 



(77) 



Since the value of the ground state energy in the 
absence of potential, i.e., for ct = 0, reduces to the 
first two terms of Eq.(77), we get the condensation 
energy for arbitrary filling, i.e., N not exactly equal 
to as 



pA^ 



(1-a) . 



(78) 



Again, this result fully agrees with the expres- 
sion of the ground state energy derived from 
Richardson-Gaudin equations for arbitrary TV and 
a, as given in Eq.(72). 

The condensation energy, given in Eq.(78), can 
be seen as each of the N correlated pairs bring- 
ing its own "binding energy" to the total conden- 
sation energy, this binding energy being linearly 
decreased compared to one isolated Cooper pair 
by a factor (1 — N/pfl) due to the Pauli exclusion 
principle between the electrons from which the cor- 
related pairs are constructed. 



VII. EXCITED STATE ENERGY 

Let us end by considering excited states in the 
case of arbitrary filling. Their energies read £b + 
i^k with i?k given by Eqs. (34,55). The lowest 
excited state energy thus corresponds to the lowest 
possible value of = (ek — p)'^ ■ 

We first note that ek is inside the energy layer 
{epoi^Fo + ^) by construction, while p, does not 
necessarily fall inside this layer. If p is inside this 
layer, it is energetically favorable for the lowest 
excited state energy to have ek = p- This leads to 
the standard result 



A. 



(79) 



The difference between the minimum excited state 
energy and the ground state energy £b is then 
equal to A, commonly called " gap" . 

In the standard half-filled configuration, ek can 
be equal to p because p = eFo +i^/2 for any a. For 
arbitrary filling however, Eq.(76) shows that there 
are two pair numbers A'^i and which are such 



that p is either below epo for N < Ni or above 
epg+n for N2 < N. For such iV's, we cannot have 
ek = These pair numbers are precisely given by 



iVi = pn 



N2 = pfl- pfl 



1 + a ' 



I+CT 



= pn-Ni 



(80) 



(81) 



Since the number of pairs possibly enjoying the 
potential is by construction smaller than pU, this 
shows that the N extensions of these two "anoma- 
lous" phases, < iV < TVi and N2 < N < pfl, 
are very small in the weak-coupling limit, pV <C 1, 
but exactly as wide. For N < Ni, the minimum 
value of i^k is reached for the ek lower boundary 
epg, while, for A'2 < N, it is reached for the upper 
boundary ek = ep„ + ft. Using Eqs. (34,75,76), this 
gives {Ek)^^^ for < iV < iVi as 



(Ek) 



(dilute) 
mill 



N 
P 



1 



(82) 



The first term in the RHS corresponds to the ki- 
netic energy increase of N pairs under a constant 
density of states p due to an excitation at ek = 
epg, while the second term is just half the binding 
energy of a single pair, as obtained by Cooper. 

Similarly, for N2 < N < pCl, the Ej^ minimum 
is given by 



(Ek) 



(dense) 
min 



N 



n- 



1 - cr 



(83) 



The first term again is the kinetic energy increase 
of A'' pairs due to an excitation at ek = e^o + 
while the second term is the same as in Eq.(82). 

The form of Eqs. (82,83) nicely shows the du- 
ality which exists between electrons and holes in 
the potential layer: bound pairs for N < Ni are 
formed out of electrons, while for N2 < N they 
seem formed out of holes. In order to excite the 
system, one has to break a pair; this explains why 
(-^^-Omin contains one half of the single pair binding 
energy (one half only because a full "pair breaking" 
means adding two excitations, not just one). The 
first terms in Eqs. (82,83) come from kinetic energy 
of the excited electrons suffering Pauli blocking. 

The existence of a duality between electrons and 
holes is also supported by the standard result of 
('^fe)min obtained for TVi < < N2, namely, 
('^fc)min = A as given in Eq.(79): the situation 
then is fully symmetrical with respect to the mu- 
tual exchange of electrons and holes. For pV <^ 1, 
i.e., in the weak-coupling limit, A^i and pU, — N2 es- 
sentially correspond to the number of states within 
half the single pair binding energy. Therefore, 
N < Ni corresponds to a "superdilute" regime 
of pairs (or a "superdense" regime of holes), with 
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wave functions of individual pairs overlapping only 
slightly, while N2 < N corresponds to a "superdi- 
lute regime of holes" (or a "superdense" regime 
of electrons). Note that, with increasing pV, the 
Ni boundary approaches prt/2 from below, while 
approaches pQ/2 symmetrically from above. 
Hence, an interval of iV's which correspond to the 
usual BCS excited state with (-Efc),!;;,! = A shrinks 
to zero when pV — >■ 00. 

We wish to stress that the transitions at = A^i 
or N2 must be smooth with respect to the ground 
state energy because the same expression (77) 
holds in the three regimes. By contrast, for the 
first excited states, these regimes correspond to dif- 
ferent kinetic energies for the excitation. As easily 
seen from Eqs.(75,82,83), this induces discontinu- 
ities in the higher-order derivatives of (-E'fc)^;^ with 
respect to pV, on both sides of the two transitions. 

The configuration considered here, with a pair 
number not exactly equal to half-filling but be- 
ing a varying parameter, can be seen as a model 
for a density-induced crossover between isolated 
fermionic molecules when N is very small, towards 
a dense regime of Cooper pairs, typical of BCS 
superconductivity. We then expect the excitation 
energy in the dilute regime of pairs to be controlled 
by the single pair binding energy, as we find, while 
at higher densities, it is controlled by cooperative 
many-body effects. For density induced BEC-BCS 
crossover, see also Ref. [21]. 

The analysis presented in this section has simi- 
larities with the BEC-BCS crossover considered by 
Eagles[18] and by Leggett[19]: Eagles keeps a non- 
constant 3D density of states, while Leggett uses 
a potential without an upper cut-off at epo + 
the irrelevant divergences being then cured by a 
"regularization" procedure. The nice aspect of the 
model we here consider, with a upper cut-off and 
a constant density of states — exact for 2D super- 
conductors — is that it leads to analytical results. 
Moreover, the upper cut-off at epo + ^ evidences 
an interesting "superdense" regime of pairs for N 
close to complete filling pfi, this regime being un- 
derstood as a dilute regime of holes in the poten- 
tial layer with similarities with the dilute regime 
of electrons when N is very small. 



VIII. CONCLUDING COMMENTS 

In this paper, we show that, through a proper 
choice of the Zk scalars appearing in the Bo- 
goliubov Hamiltonian, the difference between the 
ground state energies of the BCS and the Bogoli- 
ubov Hamiltonians can be reduced to underex- 
tensive terms, negligible in the thermodynamical 
limit. The energy of this ground state fully agrees 
with the value obtained by analytically solving 
Richardon-Gaudin equations which give the exact 



eigenstates of the BCS Hamiltonian in the canon- 
ical ensemble. Actually, this agreement is rather 
puzzling because the ground state of the Bogoli- 
ubov Hamiltonian corresponds to the wave func- 
tion ansatz proposed by Bardeen, Cooper, Schricf- 
fer, with all pairs condensed into the same state, 
while the exact Richardon-Gaudin eigenstate reads 
as a product of correlated pairs which are by con- 
struction all different due to the Pauli exclusion 
principle. 

At the present time, we have been able to get 
the Ri parameters, solutions of Richardson-Gaudin 
equations, analytically in the dilute limit only on 
the single Cooper pair scale [6, 11, 12]. These pa- 
rameters, which rule the A'^-pair wave function, as 
seen from Eqs.(l) and (2), have imaginary parts 
which run away from the real axis when N in- 
creases. Unfortunately, we are not yet able to write 
the analytical form of each Ri in the dense limit. 
Instead, we have found a way to get the sum of 
these Ri^s in the dense limit, the singular imag- 
inary parts disappearing from the sum which is 
real by construction since it gives the iV-pair en- 
ergy. The extensive term of the Ri sum we find, 
exactly matches the extensive part of the BCS 
ground state energy, as obtained through the Bo- 
goliubov approach [sec Eqs.(72) and (78)]. 

So, at the present time, we have two totally dif- 
ferent forms for the ground state wave function 
which exactly give the same extensive part of the 
ground state energy: the BCS ansatz which also is 
the Bogoliubov Hamiltonian ground state, its pro- 
jection on the A^-pair subspace reading as 

\^BCs)=B^''\0) , (84) 

where B'^ = SkC^k/wk) with i?^ = ^^k^-k' ^^^^ 
the exact Richardson-Gaudin wave function which 
reads as 

IV(iV)) = B\R,)B^{R2) ■ ■ ■ B\Rn)\0) , (85) 

where B^{R,) = Ek 5k['-'k/(2ek - R^)], the R^s 
being by construction all different. Due to Pauli 
blocking, it is clear that, when N increases, each 
of these states get more and more different from 
a bare juxtaposition of either B^ or B''{R) cor- 
related states. As a result, we cannot totaly ex- 
clude that, when TV gets very large as in the ther- 
modynamical limit, these two states end by corre- 
sponding to the same linear combination of differ- 
ent B^_^B^^ ■ ■ ■ i^kjv products. Since a precise mi- 
croscopic understanding of superconductivity goes 
through the knowledge of the correct ground state, 
to study the link between the BCS ansatz and the 
exact Richardson-Gaudin wave function appears to 
us as highly desirable. We hope that the present 
work will help to reopen a field commonly consid- 
ered as fully understood. 
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